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ABSTRACT 

Variations on the traditional cost minimization of continuous review 
formulation are investigated in an effort to improve service as measured 
in terms of time-weighted shortages per unit time. It is proposed that 
the minimization of time-weighted shortages per unit time will improve 
service in current Navy Supply Operations. Various models are presented, 
without reliance upon unknown parameters such as order cost and carrying 
cost, with necessary conditions and solution algorithms. 
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I. INTRODUCTION 



Inventories exist to provide service to customers by satisfying their 
demands from on- hand material. It follows, then, that a reasonable 
objective of inventory management is the maximization of service provided 
which is achieved by minimizing stockouts. In particular, total time- 
weighted shortages is thought to be the desired objective. 

In pursuing this objective, the manager of a realistically large, 
multi-item inventory system has a number of constraints imposed on his 
,! when to buy and how much to buy” decisions. The stock points of the 
Navy Supply System have investment and reorder workload constraints which 
are real and binding. 

The classic variable cost minimization formulation is the most used 
method for solving this inventory problem. Multi-item problems are 
usually solved by assuming that they can be dealt with as a series of 
independent single item problem. In the presence of binding constraints 
on a population of items this approach is not applicable. Additionally 
the cost minimization formulation requires the estimation of cost para- 
meters which are arbitrary or at least very difficult to estimate. 

As a consequence of this argument, a series of models are formulated 
for multi-item policies subject to investment and reorder workload constraints. 
These models do not employ the standard ordering, shortage costs. This 
approach was suggested by A. P. Tully QQ . 

In the next section, the problem formulation and the general models 
are developed. Section III develops the single item model as preparatory 
to studying multi-item case. Section IV presents a simplified multi-item 
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formulation in whicn only the items reorder points are decision variables. 
The general multi-item continuous review model is developed in Section V. 
In the final section, we summarize and state some tentative conclusions. 
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II. FORMULATION 



It is desired to formulate inventory decision rules for multi- item 
inventories subject to specific constraints. The inventory decision rules 
will be of the reorder point - reorder quantity, continuous review type. 

As suggested by the introduction, the formulation to be used involves 
the minimization of total time-weighted shortages subject to: 

(1) Total average investments costs investment limit; and 

(2) total number of orders £ reorder workload constraints. 

Note that such a formulation would not be based on minimizing variable 
costs . 

The specific form of the model depends upon the assumption about the item 
demand characteristics and expressions for the total average on-hand inventory 
level and total number of buys per unit time. The first assumption is the 
distribution of lead time demand is normal for all items. 

The following notation is used throughout the paper. For the i-th item let 
Ci = item unit cost in dollars; 

4A = mean demand per unit time in units; 

Mi = mean lead time demand in units; 

St = Standard deviation of lead time demand in units; 

$(r^) = probability that lead time demand exceeds r; 
r^ s reorder point; and 
Qi = reorder quantity. 

Also let; 

= investment limit; and 
K 2 - reorder workload constraint. 
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With a continuous review inventory policy an order is placed after 
the demand of Q units. It follows then that the average number of orders 



per unit time is . For a multi-item inventory with N items, the total 

expected number of orders placed per unit time is 




(II. 1) 



Total inventory investment is the priced-out value of the total 
expected on-hand inventory. As shown by Hadley and Whitin for 
continuous review the expected on-hand quantity, E(OH) is given by 



E(OH) = r + | -/X + B (Q,r), 

where B(Q, r) is the expression for the expected shortages at any point 
of time. If lead time demand is normally distributed it can be shown 

DO that 

B(Q, r) = ~ Cp(r) - p (r + Q) ] , (II. 2) 

where p(r) = ~ ]j 2 + (r ~/l ) 2 U $ ( 

. 1 - -Ly^ 

q>(X) = — e 2 ; and 

f2JT 

$(r) = J <|>(X) dX. 

r 

The expected on-hand quantity expression can be simplified by omitting 
the B(Q,r) term, and this approximation is reasonable if the risk of stock 
out is not too large. This assumption is employed throughout the thesis. 
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With this assumption, the total inventory investment is then given by the 
expression 



Ci ( ri + ^ 
i=l 



At). 



(II. 3) 



The expected number of backorders at any time may be explicitly 
determined from the steady state probability distribution for negative 
net inventory levels. Hadley and Whitin []2^| used this approach and 
showed that when lead time demand is normal, the time-weighted shortages 
expression is given by Eq (II. 2). 

If the risk of stockout is small, then the expression for the time- 
weighted shortages can be simplified by ignoring the |5(r+Q) term, which 
yields expected time-weighted units short per unit time for the ith item 
as : 

p(ri) (II. 4) 

The objective can now be stated as the minimization of the time- 
weighted shortages for the entire inventory, and the formulation is: 



N 

” in Z ■ p(ri) (II. 5) 

i=l 

subject to: 

N 

^ Ci(ri+^-jUi) L \ > (II. 6) 

i=l 

~ k 2 » - (ii. 7) 

i=l 

unrestricted. 



Qi y 0, and r^ 
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Note that the investment constraint will always be active, given the 
objective function used, but that reorder constraint may or may not be 
active in a given problem. 
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III. SINGLE ITEM B(Q.r) MODEL 



The basic continuous review formulation for one item with normal 



OUL, tf 2 ) lead time 


demand can be written 


as 




min: 


nr 




(m.i) 


subject to: 


C (r + | - /t ) 


4 K i » 


(III. 2) 






4 k 2 , 


(III. 3) 



Q 

Q > 0, and r unrestricted. 



A. NECESSARY CONDITIONS 

To solve the single item continuous review model, set up the 
Lagrangean function 



L(Q, r , /Jl 



e> = i 



£(r) 



+ <nC c ( r+ 2 - +0H^- K 2l] 



Taking the partial derivatives 



3L 3L QL , QL 

QQ ’ -Dr * ‘dfl ’ and *30 5 



and setting the resulting expressions equal to zero yield; 



*£L 




( r ) + 

2 Q 2 



0 
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*b L 



5 C(r-u) § <Ff) - * <(> (-^-) 3 +41C-0, 



(111. 5) 

(111. 6) 



'dL 



1L 

t»e 



C( r + | -/t ) 



K = 0 and 
1 




(111. 7) 

(111. 8) 
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These equations may be rewritten as 



Q 2 = 2 Cp(r) + M3 , CUI. 9) 

cf <rjjr) - (r-A)^ <r£) =U C Q> (m.io) 

C (r+| -A) - K lf (III. 11) 

and = K 2 . (III. 12) 



These are the necessary conditions for solution. 



B. ITERATIVE SCHEME 

If the reorder constraint is active, the reorder quantity is determined 
from the equation (HE. 12) 



Q £ 

q = A 

k 2 



which implies: 



(III. 13) 



Note that in this case, equation (III. 11) can be solved for r yielding 



r = 



2Kr 



+ A • 



(III. 14) 



Hence (Q*, r*) are uniquely determined from equations (III. 13) and (III. 14), 
If the reorder constraint is not active, the reorder quantity is 
determined along the line which is equation (III. 11). It is observed that 

j^l 

the boundary conditions are as follows: when r = 0, Q = 2 ( — + ytt ) , and 

C 

K-, 

when Q = 0, r = — + A Plotting these conditions, Figure 1 is obtained. 

C 
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is greater than 2(-jr + /4 ), then the 



It is realized that if 



K, 



investment and reorder constraints are mutually inconsistent (infeasible) 
If all solutions lie on the line C(r + ^ - JX. ) = K for Q > ^ 






and 



, Kn /i 

C 2K , then pick up a Q, solve for from the equation (III. 11), 



evaluate 



H il 



and continue to iterate on Q until 



Mil 



minimum is 



found. Here parameterizing will give shortages as a function of the 
investment limit. 

An alternative approach would be to use a Lagrangean function with the 
boundary conditions jj—* ^ Q ^2(^i ) ~| and the equation (III. 9) 

with 0=0. From these equations, it is determined that 



LX*) 

Ki 2 

2( — +/L) C 



(III. 15) 



Hence a double iteration is required as the function of /J^, Q, and r ; i.e 
iterate Q, r for a given value of the multiplier. The solution is the 
minimum time-weighted shortages for the multiplier used; minimum for the 
investment level (K^) imputed to that multiplier value ^ . Then a sug- 
gested procedure would be a binary search. 

Plotting shortages vs. multipliers on the Q, r and shortages planes, 
the different level of shortages will be obtained. 



C. EXAMPLE OF THE SINGLE ITEM MODEL 

Consider an item with its distribution of lead time demand normal 

Let P>= 125, d 2 = 25, % = 500, C = 12, = 720, and K 2 - 8. 
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The problem is 



minimize 



subject to: 



z» i p (r) 

C(r + | -/L) £ K , 
2 1 



% 

Q 



4 



K 



2 ’ 



Q>0, and r unrestricted. 

From ^ £ K^, it is determined that Q '%■ = 62.5. Hence if the 

reorder constraint is active, Q* = 62.5, r is obtained from (III. 11) as 
r * = 144, and Z* = .606. 

If the reorder constraint is not active, then a double iteration is 
required as the function/]^ , Q, and r . The solution is the minimum time- 
weighted shortages for the multiplier used. Iterating on Q, and ignoring 
the reorder constraint, produces the following results: 



Q 


r 


Time-weighted 
shortages per 
unit time 


Expected number 
of units short 
per unit time 


10 


170 


0.3198 


17.8446 


18 


166 


0.027 


14.5262 


36 


157 


0.333 


16.2599 


50 


150 


0.471 


20.8288 


62 


144 


0.606 


25.0432 


76 


137 


0.904 


32.9660 


90 


130 


1.247 


42.6242 


105 


122 


1.741 


53.6828 



Table 1. Time-weighted Shortages and Expected Number of Units 
Short for Single Item Case. 

Note: Two approaches seem to follow together. 
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IV. SIMPLIFIED MULTI- ITEM B(Q,r ) MODEL 



The basic continuous review formulation for the multi-item problem 
was given as: 



minimize Z = 



I 

i=l 



subject to: N 



Qi 



Ci( ri + ^ - At) 4 K x , 



(IV. 1) 



(IV. 2) 



i=l 



N 

I 

i=l 



— 4k 

Qi = K 2 ’ 



(IV. 3) 



Qi > 0 , and r^ unrestricted. 

Suppose the order quantities are fixed by some other criterion. 
Specifically the assumption is made that order quantities are determined 



VU 

from the equation Qi = G I — > 



(IV. 4) 



where G is a constant, which is assumed to be the same for all items. 

From the reorder constraint (IV. 3) and the assumption (IV. 4) it follows 

that 



or 




(iv. 5) 



The determination of G then fixes the order quantities from equation (IV. 4) 
and eliminates one set of decision variables from the problem, i.e.. 
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Qi = 



i=i 



K„ 




(IV. 6) 



Substituting (IV. 4) into equation (IV. 2), it follows that 



N 



E 

i=i 



C1 < ri + T j¥ ■ jLv > ^ \ 



Equation (IV. 7) can be reduced to the form 



(IV. 7) 



N N N 

qr^^. y y jfoCi + £cij*i = . 



(iv. 8) 



1=1 i=l 1=1 

Now the multi-item B(Q, r) Model with fixed Q's can be written as 



N 



minimize Z = 



subject to: 



E 

i=l 



ft(ri) 

Qi 



N 

Z c i r i 4 k ; . 

i=l 1 



and ri unrestricted, 

N 

where Qi = 



2K„ 






y Ci 



= G 



EH 

\jcT~ 



A. NECESSARY CONDITIONS 

To solve the simplified model, set up Lagrangean function 



L ( 



ri ’^ " f “g" fi( rl) Ct rt " K x j ( IV - 9 > 
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Taking the partial derivatives with respect to decision variables and 
setting these expressions equal to zero yields 



j5l_ _ JL (cT 
0ri G Jfl.i 



[(ri-A)$(- ^~ ) - *%Ct = 0, 



and 



2^ = 




Ci ri 




Thus the conditions for solution to the problem are 




and 




Ci ri = K ’ 



(IV. 11) 



B. NUMERICAL SOLUTION SCHEME 



In general, these equations cannot be solved in closed form. A numerical 
solution procedure is suggested. Let us consider a numerical method of 
solving equations (IV. 10) and (IV. 11). Observing equation (IV. 10), which 
is 




1 

gJ Ci'Xi 



[<A <h- £ w L >- c 




ri-Ai 
~an — 




note that the right-hand side has a lower bound of 0 since %>, 0 from 
G > 0, Ci£i>0 and 0*i <j> ( ^— ) - (ri-Ai) $ ( 0 for all ri. 



Specifically^ = 0 when 



r i becomes infinite and 






§ (^) - o, bo. 



r^ =00 violates the investment constraint. 
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Hence for the initial value of^ it is reasonable to start with 



r^ = 0, then this implies 



% ~ T== Nit > + 

g L 



jUi 

<fi 






for all i, 



or 



*' v '“" j= 



Let 



for all i ; 



= § implies that there is at least one r at its lower bound of zero. 

C 

Hence - — is a convenient starting point. Then a suggested solution 

2 0 

procedure would be to begin at /J^ = — , solve equation (IV. 10) for the 

r^'s and compute the value of constraint using equation (IV. 12), which is 



N 



E 

i=l 



C i r i "H. 



(IV. 12) 



* b 

A bisection search will be used again. If K , increase^ by — — . 
If H ^ K^ f , decrease ^ by -iL- , Recompute the r^'s and the value of 
constraint using equation (IV. 12). If the increase (or decrease) of 
has not caused the change of inequality sign, increase (or decrease) 



by the same amount 



If the sign of inequality has changed, then 



reduce the increment to and increase (or decrease)^ , solving for 
the r^ f s at each value of 4^ and computing the value of H until the 
sign of the inequality switches again. Continue until H = K^ 1 or until 
H is within some tolerable limit of K^ f . This approach will converge to 
the optimal solution rapidly. 
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From the Kuhn- Tucker theorem {j3T\ , if we have a convex objective 
function and a convex constraint region, the necessary conditions are also 
sufficient. Since the constraint under consideration is linear in r , the 
region is convex. To show Z( r i) is convex, consider the equation of the 
expected time-weighted shortages, 



Zi - 1 



Ci ^(ri). 



G 

2 

Now if ^-5 — ^ 0 for all ri, then Zi i s convex. Taking partial 

'brt 2 

derivatives , 



DZi 



and 



*0ri 



% 

QZ i 

< dr i 2 



<rj| C ( ^ °- 



_1_ 

G 



c i 5<-^) • 






(IV. 13) 



Equation (IV. 13) will always be greater than or equal to zero. Under 
these conditions Zi is convex. It follows that Z is convex since it is 
the sum of convex functions. 

C. EXAMPLE OF THE SIMPLIFIED MULTI- ITEM B(Q,r ) MODEL 

Let us consider an inventory of three items. It is assumed that the 
distribution of lead time demand is normal with mean A i and variance ^ i 2 
for the i*-* 1 item. Let the item data be as follows: 





Item 1 


Item 2 


Item 3 


Ok 


1000 


1500 


2000 


Ci 


1 


10 


20 




100 


200 


300 


O'i 2 


100 


100 


200 



K 



l = $8,000 and K 2 = 15. 
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Now the problem is 



minimize Z 



3 




subject to: 



Ci ri 4 = K x - -y 

i=l 



3 

y~ j^ici + 

i=l 




Ci Jit , 



where 

G 




= 23.6065. 



order quantities are determined from equation (IV. 4) , 



Q1 = 746.5022, 

Q2 = 289.1190, and 
Q3 = 236.0648 . 

With these order quantities, the problem now can be stated as 



minimize Z 



3 

4 =1 



1 

23.6065 




p ( ri) 



subject to: 




Ci ri 4 8000 + 



3 

Z 

4=1 




11,920.5078. 



Notice that the solution to the above problem is that vector r such that 

4 = —1 |Vi <j>(— j yffi ) - ( ri-Ai)$ ( - Vf'^ 1 )] » for a11 (IV. 10 ) 

G Jclri L J 
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For each value of , there will be some (. r } for which equation 
(IV. 10) is satisfied. However, having the convex objective function, 
there exists only one £ri} for all i such that equations (IV. 10) and 

g 

(IV. 11) are satisfied at the same time. For the initial value of ^ 9 

using 8 = — min „ * [V;<^(- 4^) + § (- — -)1 

G L di <ri J 

for all i when ri = 0, ^lo = 0.0318 will be used. Decreasing from 
0.0318, 'ftp = 0.0055 is obtained, which implies 

r l = 234.3750, 
r 2 = 264.0625, and 
r 3 = 453.1250. 

Checking the constraint, it is found that 



^ Ci ri = 11937.5000, 

1=1 

which is within 0.2 percent level of k| 



is 



The expression for the time-weighted shortage per item per unit time 

*1 'ijf [H°’ l 2 +<ri ' A) 2 i^ < ' £ ^ a) 'i • 

The time-weighted shortage per unit time is computed as 

Z 1 ■ 23^065 [ j { 100 2 + T 2 <«°- l 00)4O] 



= 0.2179 



Z 0 = — 
2 23 



nrr.r 

.6065 J 1500 [ 2 l 



100 2 + (285-200) (285-200) 



<K §>] 



= 2.7634 , and 
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are 



Zq = 



23.6065 J2000 



J2000 l 2 [ 



200 2 + (285-300) 2 



}$<- 



15 

200 



) - 



200 



= 10.5232 



(285-300) 




Total time-weighted shortages per unit time for the entire inventory 



3 

Z = ^ Zi = 13.5045. 
i=l 
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V. GENERAL MULTI- ITEM B (Q , r ) MODEL 



Let us consider the completely general continuous review model for 
multi-item B(Q, r ) Model which was: 

N 



minimize 



subject to: 



Z = 



N 



E 

1=1 



-Si- p(ri) , 
Q)- 

Si 



Ci ( ri+ 2 - P* i) 4: K 
i=l 1 



N 

E 

i=l 



'k. 

Qi 



£ K 



Qi y 0, and unrestricted, 



where 

and 



f 5 ( ri) - I [°i 2 + ( ri-Ai) Z J 5 ( --■ £- ) - % ( ■ 



Wi is a weighting factor, 



A. NECESSARY CONDITIONS 

Applying the Lagrange Multiplier technique, the following results will 
be obtained: 



L(Qi, *i>%> ®) =' 



N 



3Qi QJ 



3 L _ (tiit 
*0ri Q 

2JL = ^- 

2_ 

i=l 






Ci 



P(ri)+,/y" Ci(ri+ S | 


1 fN 

-Jiiy d « Vqi-v (v.d 


h=i 


J l i=1 3 


Cl? 

II 

o 


(V . 2) 




-HjJCi = 0; (V . 3) 


Ai) - K =0; and 


(V.4) 
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0. 



(V . 5) 



3l_ = VAi- 
29 l— Qi 2 

i=l 



Then equations (V.l) and (V.2) can be solved for each multiplier yielding 



6 = 


- ( ri-/ii)$ i ~) - 2W^( ri) 






2 'Xi W i 




11 

CM 

O' 

u 

o 


2 fwi £( ri) + 9%i] 
%Ci 


(V.6) 


and 

% - 


«i$( "Vi" 1 ) - < « 


(V .7) 


CiWiQi » 


N 

E 

i=l 


Qi „ 

Ci( ri + 2 -A i) = K , 


(V.4) 



and 



N 



I 

< =1 





(V.5) 



Considering equation (V.7), the right-hand side of equation is always 
positive, which implies 0, since ^=0 implies Qi = CO which 

violates equations (V.5) and (V.4) . 

Observing equation (V.6), the right-hand side of this equation is 
always positive. This suggests two possible cases for considerations. 

Case I . %>0, 0- 0 (Ignoring the reorder workload constraint). 

The necessary conditions in this case are: 

2 _ 2 Wi P( rl) 

Q i /^Ci i ( V * 8 ) 
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(V .7) 




dt <(> (_£^«i) - ( ri -Ai) $ ( ffi 1 ) 

CiWiQi 



> 



and 



Ci ( r i + ^ " Ai) 
i=l 




(V.4) 



Case II , %> o, 0 > 0 (Both constraints are active) . 

The necessary conditions are the same as the previous page with Eq 
(V.6), Eq (V.7), Eq (V.4), and Eq (V.5). 



B. ITERATIVE SCHEME FOR CASE I 

Equations (V.8) and (V.7) can be thought of describing two curves in 
the Q-r plane. Equation (V.8) clearly shows that order quantity approaches 
zero with increasing reorder point to infinity. Furthermore, from equation 
4^ ^ 0 for all i; 

- [ <^> - < rl- 



uv j. 

(V. 8 ), drT 



m. = 

dr,- 



Taking the second partial. 



J 2 4^CiWip( ri) 

is greater than zero or less than zero 



d 2 Qi 

dr. 2 



depending upon the sign of numerator, which is 

d 2 Ql = 2fyCiWi p(ri) $ ( r ^~ < ^ i ) - - (ri- Jl^ (^-1)] 

pz^ci Wi J3 (ri) J \ 



dr i 2 



From equation (V.7), note that the nonnegative quantity 



<{> ( ^£i) . (rl .M) $ ] 

decreases to zero as r^ increases to infinity, at which time the order 
quantity will be zero. However, taking the first and second partials 
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from equation (V.7), 



dri 

dQi 



^ 0 will be obtained; 



< 0. 



aq 



ri 

2 



drl = . ^CiWl_ / 0 for all t 

dQl 4<^> 



n 

and d ri = 4LciWi O' i y 0 

aqi 2 " <t‘< £j $£ L ) 



for all i . 



If one plots the two curves, Figure 2 will be obtained. 

By setting the value of ri to zero in Eq (V.7) and Eq (V.8), Qa and 



Qg will be obtained respectively, which are 

- (ri - Ai) $ 

q a - — — 

% CiWi 

and 

1 2Wl[j( </l 2 +Ai) 2 $<- ^f-) + ^A<}> (- -^f)] 
Qb = -Ilci 



To have a solution, it is necessary to have Equating Eq (V.7) 

and Eq (V.8), yields 



Qi 



2 



2Wi ft(ri) 

^Ci 



[ o'i 4 ( . r -^f‘ ) - 

(^CiWi)2 



or 




2CiWi 3 p(r.) 



(V.9) 



Notice that Eq (V.9) is the function of the decision variable r only. A 
double iteration as a function oftf) , Q and r is required, i.e, iterate 
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Figure 2. Q as Function of r for Case I 
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Q, r for a given value of the multiplier . For the initial value of , 
ut 



8 = min 



2C i W i (r.=0) 



for all i, pick 



as a convenient starting point. Then solve 



equation (V.9) for the £r^’s} , compute £q^*sJ using Eq (V.8), and 

N Q • 

calculate the value of constraint using Eq (V.4). Let 



J2 c i( r i + -r-Ai>=H. 

i=l 



A binary search will be used. If H^K^, increase by 



If H < K , 



decrease 41 by 



Compute the value of H. If the increase (or decrease) 



of /]^ has not caused the sense of inequality sign, increase (or decrease) 

g 

by the same amount — . If the sign of the inequality has changed, then 

C 

reduce the increment to , and increase (or decrease)^ , solving for the 

8 

£r^ ! s J and £ Q i *s}at each value of ^ and computing the value of H until 

the sense of inequality changes again. Continue until H = or until H is 
within the tolerable region of . This method will converge to the solution 
vector rapidly. 

From the Kuhn- Tucker theorem QT] , if we have a convex objective 
function and a convex constraint region, the necessary conditions are also 
sufficient. Since the first constraint under consideration is linear in 
r, the region is convex. To show Z(r^) is convex, it is required that 
the Hessian is greater than zero or equal to zero; i.e., 

.2, 



2_ 

QZi 






'dr i 'dr. 



2 

•dZi 



< 3Qi < 3 r i *dQid Qj 

Let oC Oi) = d (-■ — ) - ( r i 



^ 0 
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Taking first and second partial with respect to each decision variable, 
the following will be obtained: 



3Zj. = _ <* <r t ) 

O r i Qi _ 



: = 



*d 2z i 

r. 



Qi 



$ 



r i-Ai 

<fi 



Qi 2 



'i z i 



*3 Z i _ 0<( r i) 



Qi 2 ’ 



^i^Qi 
From Equation (V.10), 
^ Z i _ 



and 



t£z t 



Qj© Qj 

= < r i> 



2 P <gj> 
Qi 3 



"dQlB r t 



t) 2 Zi 



3 2 Z = $ ( flTj ) 2 ^ r j) . <=* < r i> °t< r j> 



rj Qj 9*^9^ *9^9 Qi Qi Qi 3 <^2 q 2 



or 






i > 

(V.ll) 



Hence Equation (V.ll) is greater than zero or less than zero depending 
upon the sign of bracketed quantity. Under these considerations, the 
convexity of Zi is unknown; only a local minimum can be assured. 



C. EXAMPLE OF THE GENERAL MODEL FOR CASE I 

Once again consider the inventory of three items from IV. C, it is 

assumed that the distribution of lead time demand is normal with mean 

2 th 

Mi and variance di for the i item. 





Item 1 


Item 2 


Item 3 


u 


1000 


1500 


2000 


c i 


1 


10 


20 


At 


100 


200 


300 


di 


100 


100 


200, and let 



K x = $8,000 and K 2 = 15. 
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Reviewing section (V.A), note that the solution vectors must satisfy 



q 2 _ 2 »! 
i fti r 



/rt r 



(V.8) 




(V. 7) 



and 




(V.4) 



i=l 



Using the search scheme described in section (V.B), this example was 
solved with the results shown in Table 1 utilizing three investment 
levels. Comparison between the General Model for Case I and Simplified 
Multi-item Model of section III was made for the purpose of determining 
how good or bad the simplified model performed in terms of time-weighted 
shortages and marginal cost imputed to the first multiplier. These results 
are shown in Tables 2 and 3. It can be seen that the Simplified Model 
produced time-weighted shortages which were from 25% to nearly 100% larger 
than the General Model for Case I. 

D. ITERATIVE SCHEME FOR THE GENERAL MODEL 

Again equation (V.6) and V.7) can be thought of as describing two 
curves in the Q-r plane. Equation (V.6) shows that the order quantity 
approaches a fixed value, Q L as the reorder point goes to infinity. 






^2-JiC.Wi [fkr^+e^J 
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ITEM 


1 


2 


3 


REMARKS 


VARIABLE 










r 


340.6250 


342.1875 


506.2500 


K x = 8,000 


Q 


60.8684 


79.4553 


178.4906 




% 


0.0044 


0.0044 


0.0044 




z i 


0.1343 


1.7529 


7,8755 


Z = 9.7627 



r 


291.4063 


264.8438 


307.8125 


K x = 4,000 


Q 


69.4359 


101.1958 


247.2386 


% 


0.0154 


0.0154 


0.0154 


Zi 


0.5334 


7,7740 


37.9866 


Z = 46.2940 



r 


400.0000 


425.0000 


687.5000 


K x = 12,000 


Q 


52.6381 


60.4982 


137.0932 


U 


0.0007 


0.0007 


0.0007 


Zi 


0.0194 


0.2229 


1.0100 


Z = 1.2522 



Table 2. The General Model With nn Inactive Reorder 
Workload Constraint (Case I) 
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ITEM 

VARIABLES. 


1 


2 


3 


REMARKS 


r 


234.3750 


264.0625 


453.1250 


K l = 8,000 


Q 


746.5022 


289.1190 


236.0648 


K 2 = 15 


41 


0.0055 


0.0055 


0.0055 




z 


0.2179 


2.7634 


10.5232 


Z = 13.5045 


i 











r 


175.7813 


181.2500 


296.8750 


K x = 4,000 


Q 


746.5022 


289.1190 


236.0648 


K 2 = 15 


41 


0.0173 


0.0173 


0.0173 




z i 


0.8319 


11.4189 


42.8919 


Z = 55.1428 



r 


300.0000 


343.7500 


612.5000 


K 1 = 12,000 


Q 


746.5022 


289.1190 


236.0648 


m 

i— i 
1! 

CM 


% 


0.0011 


0.0011 


0.0011 




h 


0.0386 


0.4632 


1.6465 


Z = 2.1483 



Table 3. Simplified Multi-item Model Solutions 
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The second partial, 



2 

d Z Qi 



dr j 



2 “ , is greater than zero or less than zero 



depending on the sign of the numerator; 

d 2 Qi 2-fyCiWi [(3(r. ) + Ofti] $ (r t ) 2 

dr i 2 |*2^CiWi [ f( ri ) + 0?-i. } J | 

Observing equation (V.7), note that the same properties as section 
(V.B)will be obtained, which are 



dT 4 < 0, and > 0 for a11 i ‘ 



dQi 



dQi' 



If one plots the two curves something like Figure 3 will be obtained. 

Observing Figure 3, the numerical search in Hadley and Whitin Ql[] 

will be used for the general model. To initiate the numerical procedure, 

it is necessary to determine reasonable multiplier values. These two 

fo/Vi 

multipliers are used to compute r when Q T = — = Q. The solution 

c 

could be started with the point where Q = Q on the curve defined by Eq 

Li 

(V.6). The r value so obtained is used in Eq (V.7) to compute a new Q 
value Q 2 ; i.e., move from the point (Q^, r-^) on the curve defined by Eq 
(V.6) to a point on the curve defined by Eq (V.7) having the ordinate r^. 
The Q 2 value is used in Eq (V.6) to compute a new r; i.e., we move from 
the curve defined by Eq (V.7) to the curve defined by Eq (V.6) at constant 
Q. Hence a series of steps is obtained as shown in Figure 3. It is 
clear that the numerical method must converge to Q* and r* for every i, 
given Q a ^ Q b . 

Then compute the value of Eq (V.5) and Eq (V.4). 
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In changing the values of two multipliers sequentially, ^ = 0.01, 

0.02, ... and 0 = 1,2, ... a set of solution vectors, { Q, r} will 
be obtained. Then compute the time-weighted shortages and the value of 
Eq (V.5) and Eq (V.4). ’Good* solution vectors can be found with small 
variance of Eq (V.5) and Eq (V.4) and approximate the minimum time- 
weighted shortages. Once obtaining this, iterate with more careful 
selection of multipliers, compute the objective function and check the 
variances of Eq (V.5) and Eq (V.4). Continue iterating on the multipliers 
until tolerable differences between the right and left sides of Eq (V.5) 
and Eq (V.4). 

More sophisticated approaches are available. Among these is a search 
technique proposed by Fiacco and McCormick , and a two state variables 

optimal allocation using dynamic programming. 

The sequential unconstrained minimization technique Q4 J , SUMT, is 

_ N 

based on the minimization of a new function P(X,y) = f(X) +p ^"* r l/g^(X) 

over a strictly monotonic decreasing sequence of values {M • 

Under certain restrictions that will be reviewed subsequently, there 



{x(ik>l 



that respectively minimize 



exists a sequence of feasible points 
{p<x , p k )} , and it follows that X ( P k)~> X , a solution of original 
function as > 0 (k-^Oo ). The following is a concise summary of 

the steps describing the computational algorithm: 

1. Select a point X° interior to the feasible region. 

2. Select f^, the initial value of P using 

Pi = - V f <X°) T v P< x °) / I VP( X °)| 2 > 
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or 



1 

2 



Pi 



Vfq 0 )^ 1 y f(x°) 
'Vp(x°) t h 2 1 Vp(x°) 






where p(X) = '/ g^(X), and H 2 are the Hessians of f(X) and 

p (X) respectively* 



3. Determine the minimum of P(X,f^) for the current value of 
using Gradient Methods, 



2 1 m 1 

X = X - 07 P (X ), 

or 

X 2 = x 1 - QJ^pCx 1 ) [flx.BXjJ y PCX 1 ) 

4. If k^{ , estimate solution using extrapolation formula. 

5. Terminate computations if final convergence criteria are satisfied, 
f(X) - c[x(p), Jl(?) ]< £ * The theoretical optimum value v Q is bounded 
by the dual and primal function values respectively, 

G0p),/l(e)]4v 0 4 fCx(f)] . 

If the value v Q is not within the bounds above, go to step (6). 

6. Select Pfc+1 = P]c|c > where C> 1. 

7. If 1, estimate minimum for reduced P - value, using an extra- 
polation formula. 



E . EXAMPLE OF THE GENERAL MODEL (CASE II) 

Consider an inventory of three items with the following characteristics 





Item 1 


Item 2 


Item 3 


^i 


1000 


1500 


2000 


Ci 


1 


10 


20 


A 


100 


200 


300 


A 2 


100 


100 


200 



K 2 = 15 



= $ 8,000 
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The search technique proposed by Fiacco and McCormick J_4_| was used. 
The following initial feasible solution was used: = 600, r^ = 200, 

Q2 = 270, r 2 - 260, Qj - 300, r 3 - 400. 21 - fo rU , 17 . 8084 and 

i=T Qi 



P^ = - 1.6991 were obtained. 

This example was solved with the following results, 



*1 

Q 2 

*3 

and the value 



= 362.718, 

= 266.403, 

= 302.6870, 



r = 322.1433, 
r 2 = 272.0930, 

r 3 = 425.7958, P = - 0.00002231, 



of the objective function (Z) was Z = 13.39781. 
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VI. SUMMARY AND CONCLUSIONS 



It is apparent that the general models proposed do not represent the 
ultimate answer in multi- item inventory theory. Since the conditions of 
concavity and convex sets are not maintained for the general model, the 
computational algorithms which are based on concavity and convex sets of 
objective and constraing respectively could not apply directly. Though 
the P function in section (V.D) appears prohibitively difficult to work 
with computationally, it has been understood that it is minimized efficiently 
and accurately for the great preponderance of problems solved to date, by 
means of the second order optimum gradient method. 

From the example problem we see the general problem in section (V.E) 
provides solution which is better than the solution of the simplified 
continuous review model in section (IV. C) , The major advantage of the 
simplified continuous review model is its computational ease. 

While an efficient algorithm for solution of the general problem has 
not been presented, the advent of high speed computer has opened this field 
of numerical iterative procedures for large inventory systems. 

Since the value of the constraints are more easily determined than the 
order cost and holding cost, the model proposed seems much more appropriate 
than the traditional variable cost minimization models. 

From the example problem given in section (III), it is apparent that 
time-weighted shortages and expected number of units short per unit time 
seem to follow each other closely. 

Finally, all the solutions to example problems were obtained by 
utilizing a digital computer. 
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